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We report the results of a joint analysis of data from BICEP2/ Keck Array and Planck. BICEP2 
and Keck Array have observed the same approximately 400 deg 2 patch of sky centered on RA 
Oh, Dec. —57.5°. The combined maps reach a depth of 57nKdeg in Stokes Q and U in a band 
centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 
to 353 GHz, but much less deeply in any given region (1.2 p K deg in Q and U at 143 GHz). We detect 
150x353 cross-correlation in B-modes at high significance. We fit the single- and cross-frequency 
power spectra at frequencies > 150 GHz to a lensed-ACDM model that includes dust and a possible 
contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio 
r), using a prior on the frequency spectral behavior of polarized dust emission from previous Planck 
analysis of other regions of the sky. We find strong evidence for dust and no statistically significant 
evidence for tensor modes. We probe various model variations and extensions, including adding a 
synchrotron component in combination with lower frequency data, and find that these make little 
difference to the r constraint. Finally we present an alternative analysis which is similar to a map- 
based cleaning of the dust contribution, and show that this gives similar constraints. The final result 
is expressed as a likelihood curve for r, and yields an upper limit ro.os < 0.12 at 95% confidence. 
Marginalizing over dust and r, lensing B-modes are detected at 7.0 cr significance. 

PACS numbers: 98.70.Vc, 04.80.Nn, 95.85.Bh, 98.80.Es 


I. INTRODUCTION 

The cosmic microwave background (CMB) [1], is an 
essential source of information about all epochs of the 
Universe. In the past several decades, characterization 
of the temperature and polarization anisotropies of the 
CMB has helped to establish the standard cosmological 
model (ACDM) and to measure its parameters to high 
precision (see for example Refs. [2, 3]). 

An extension to the standard big bang model, inflation, 
postulates a short period of exponential expansion in the 
very early Universe, naturally setting the initial condi¬ 
tions required by ACDM, as well as solving a number of 
additional problems in standard cosmology. Inflation’s 
basic predictions regarding the Universes large-scale ge¬ 
ometry and structure have been borne out by cosmo¬ 
logical measurements to date (see Ref. [4] for a review). 
Inflation makes an additional prediction, the existence 
of a background of gravitational waves, or tensor mode 
perturbations [5-8]. At the recombination epoch, the in¬ 
flationary gravitational waves (IGW) contribute to the 
anisotropy of the CMB in both total intensity and linear 
polarization. The amplitude of tensors is conventionally 
parameterized by r, the tensor-to-scalar ratio at a fidu¬ 
cial scale. Theoretical predictions of the value of r cover 
a very wide range. Conversely, a measurement of r can 
discriminate between models of inflation. 

Tensor modes produce a small increment in the tem¬ 
perature anisotropy power spectrum over the standard 
ACDM scalar perturbations at multipoles l < 60; mea¬ 


suring this increment requires the large sky coverage tra¬ 
ditionally achieved by space-based experiments, and an 
understanding of the other cosmological parameters. The 
effects of tensor perturbations on B-mode polarization is 
less ambiguous than on temperature or .E-mode polar¬ 
ization over the range £ < 150. The E-mode polarization 
signal produced by scalar perturbations is very small and 
is dominated by the weak lensing of E-mode polarization 
on small angular scales, making the detection of an IGW 
contribution possible [9-12], 

Planck [13] was the third generation CMB space mis¬ 
sion, which mapped the full sky in polarization in seven 
bands centered at frequencies from 30 GHz to 353 GHz 
to a resolution of 33 to 5 arcminutes [14, 15]. The Planck 
collaboration has published the best limit to date on ten¬ 
sor modes using CMB data alone [3]: r 0 .oo 2 < 0.11 (at 
95% confidence) using a combination of Planck, SPT and 
ACT temperature data, plus WMAP polarization, al¬ 
though the Planck r limit is model-dependent, with run¬ 
ning of the scalar spectral index or additional relativistic 
degrees of freedom being well-known degeneracies which 
allow larger values of r. 

Interstellar dust grains produce thermal emission, the 
brightness of which increases rapidly from the 100- 
150 GHz frequencies favored for CMB observations, be¬ 
coming dominant at > 350 GHz even at high galactic lat¬ 
itude. The dust grains align with the Galactic magnetic 
field to produce emission with a degree of linear polariza¬ 
tion [16]. The observed degree of polarization depends on 
the structure of the Galactic magnetic field along the line 
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of sight, as well as the properties of the dust grains (see 
for example Refs. [17, 18]). This polarized dust emission 
results in both E-mode and B-mode, and acts as a poten¬ 
tial contaminant to a measurement of r. Galactic dust 
polarization was detected by Archeops [19] at 353 GHz 
and by WMAP [2, 20] at 90 GHz. 

BICEP2 was a specialized, low angular resolution ex¬ 
periment, which operated from the South Pole from 2010 
to 2012, concentrating 150 GHz sensitivity comparable 
to Planck on a roughly 1 % patch of sky at high Galac¬ 
tic latitude [21]. The BICEP2 Collaboration published a 
highly significant detection of B-mode polarization in ex¬ 
cess of the r=0 lensed-ACDM expectation over the range 
30 < l < 150 in Ref. [22, hereafter BK-I], Modest evi¬ 
dence against a thermal Galactic dust component dom¬ 
inating the observed signal was presented based on the 
cross-spectrum against 100 GHz maps from the previous 
BICEP1 experiment. The detected B-mode level was 
higher than that projected by several existing dust mod¬ 
els [23, 24] although these did not claim any high degree 
of reliability. 

The Planck survey released information on the struc¬ 
ture of the dust polarization sky at intermediate lati¬ 
tudes [25], and the frequency dependence of the polar¬ 
ized dust emission at frequencies relevant to CMB stud¬ 
ies [26]. Other papers argued that the BICEP2 region 
is significantly contaminated by dust [27, 28]. Finally 
Planck released information on dust polarization at high 
latitude [29, hereafter PIP-XXX], and in particular ex¬ 
amined a field centered on the BICEP2 region (but some¬ 
what larger than it) finding a level of polarized dust emis¬ 
sion at 353 GHz sufficient to explain the 150 GHz excess 
observed by BICEP2, although with relatively low signal- 
to-noise. 

Keck Array is a system of BICEP2-like receivers also 
located at the South Pole. During the 2012 and 2013 
seasons Keck Array observed the same field as BICEP2 
in the same 150 GHz frequency band [30, hereafter BK- 
V]. Combining the BICEP2 and Keck Array maps yields 
Q and U maps with rms noise of 57 nK in nominal 1 deg 2 
pixels—by far the deepest made to date. 

In this paper, we take cross-spectra between the joint 
BICEP2/AecA,' maps and all the polarized bands of 
Planck. The structure is as follows. In Sec. II we describe 
the preparation of the input maps, the expectations for 
dust, and the power spectrum results. In Sec. Ill the 
main multi-frequency cross-spectrum likelihood method 
is introduced and applied to the data, and a number of 
variations from the selected fiducial analysis are explored. 
Sec. IV describes validation tests using simulations as 
well as an alternate likelihood. In Sec. V we investi¬ 
gate whether there could be decorrelation between the 
Planck and BICEP2 /Keck maps due to the astrophysics 
of dust and/or instrumental effects. Finally we conclude 
in Sec. VI. 


II. MAPS TO POWER SPECTRA 


A. Maps and preparation 


We primarily use the BICEP2/Aecfc combined maps, 
as described in BK-V. We also use the BICEP2-only 
and Keck- only maps as a cross check. The Planck maps 
used for cross-correlation with BICEP2/Aecfc are the full- 
mission polarized maps from the PR2 Planck science re¬ 
lease [31] [32], a subset of which was presented in PIP- 
XXX. We compute Planck single-frequency spectra as the 
cross-power spectra of two data-split maps, in which the 
data are split into two subsets with independent noise. 
We consider three data split maps: (i) detector-set maps, 
where the detectors at a given frequency are divided into 
two groups, (ii) yearly maps, where the data from the 
first and second years of observations are used for the 
two maps, and (iii) half-ring maps, where the data from 
each pointing period is divided in halves. To evaluate un¬ 
certainties due to Planck instrumental noise, we use 500 
noise simulations of each map; these are the standard set 
of time-ordered data noise simulations projected into sky 
maps (the FFP8 simulations defined in Ref. [33]). 

While the Planck maps are filtered only by the instru¬ 
ment beam (the effective beam defined in Refs. [34] and 
[35]), the BICEP2/Aecfcmaps are in addition filtered due 
to the observation strategy and analysis process. In par¬ 
ticular, large angular scales are suppressed anisotropi- 
cally in the BICEP2/ Keck mapmaking process to avoid 
atmospheric and ground-fixed contamination; this sup¬ 
pression is corrected in the power spectrum estimate. 
In order to facilitate comparison, we therefore prepare 
“Planck as seen by BICEP2 /Keck' maps. In the first 
step we use the ANAFAST, ALTERALM and SYNFAST rou¬ 
tines from the HEALPix [36] package [37] to resmooth the 
Planck maps with the BICEP2/Aec/c beam profile, as¬ 
suming azimuthal symmetry of the beam. The coordi¬ 
nate rotation from Galactic to celestial coordinates of 
the T, Q, and U maps is performed using the ALTERALM 
routine in the HEALPix package. The sign of the Stokes 
U map is flipped to convert from the HEALPix to the 
IAU polarization convention. Next we pass these through 
the “observing” matrix R, described in Section VI.B 
of BK-I, to produce maps that include the filtering of 
modes occuring in the data processing pipeline (includ¬ 
ing polynominal filtering and scan-synchronous template 
removal, plus deprojection of beam systematics). 

Figure 1 shows the resmoothed Planck 353 GHz T, Q, 
and U maps before and after filtering. In both cases 
the BICEP2/Aecfc inverse variance apodization mask has 
been applied. This figure emphasizes the need to ac¬ 
count for the filtering before any comparison of maps is 
attempted, either qualitative or quantitative. 
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FIG. 1. Planck 353 GHz T, Q, and U maps before (left) and after (right) the application of BICEP2 /Keck filtering. In 
both cases the maps have been multiplied by the BICEP2 /Keck apodization mask. The Planck maps are presmoothed to 
the BICEP2 /Keck beam profile and have the mean value subtracted. The filtering, in particular the third order polynominal 
subtraction to suppress atmospheric pickup, removes large-angular scale signal along the BICEP2 /Keck scanning direction 
(parallel to the right ascension direction in the maps here). 


B. Expected spatial and frequency spectra of dust 

Before examining the power spectra it is useful to re¬ 
view expectations for the spatial and frequency spectra 
of dust. Figure 2 of PIP-XXX shows that the dust BB 
(and EE) angular power spectra are well fit by a sim¬ 
ple power law Vi oc f -0 42 , where T>t> = C(.l (£ + 1) /2i r, 
when averaging over large regions of sky outside of the 
Galactic plane. Section 5.2 of the same paper states that 
there is no evidence for departure from this behavior for 
1% sky patches, although the signal-to-noise ratio is low 
for some regions. Presumably we expect greater fluctua¬ 
tion from the mean behavior than would be expected for 
a Gaussian random field. 

The spectral energy distribution (SED) of dust polar¬ 
ization was measured in Ref. [26] for 400 patches with 
10° radius at intermediate Galactic latitudes. Figure 10 
of this reference shows emprically that the mean polar¬ 
ized dust SED is described by a simple modified black- 
body spectrum with T<j = 19.6 K and /?d = 1.59 ±0.17 to 


within an accuracy of a few percent over the frequency 
range 100-353 GHz. Within this frequency range varia¬ 
tions in the two parameters are highly degenerate and 
the choice is made to hold fixed at the value obtained 
from a fit to the SED of total intensity, and describe any 
variation with the /3d parameter. The uncertainty on /3d 
quoted above is the ler dispersion of the individual patch 
measurements, and hence is an upper limit, since some 
of the fluctuation is due to noise rather than real vari¬ 
ation on the sky. This SED is confirmed to be a good 
match to data when averaging over 24% of the cleanest 
high latitude sky in Fig. 6 of PIP-XXX. 


C. Power spectrum estimation and results 

The power spectrum estimation proceeds exactly as in 
BK-I, including the matrix based purification operation 
to prevent E to B mixing. Figure 2 shows the results 
for BICEP2/Keck and Planck 353 GHz for TT, TE, EE, 
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and EE. In all cases the error bars are the standard 
deviations of lensed-ACDM+noise simulations [38] and 
hence contain no sample variance on any other compo¬ 
nent. The results in the left column are auto-spectra, 
identical to those given in BK-I and BK-V—these spec¬ 
tra are consistent with lensed-ACDM+noise except for 
the excess in BB for £ < 200. 

The right column of Fig. 2 shows cross-spectra between 
two halves of the Planck 353 GHz data set, with three dif¬ 
ferent splits shown. The Planck collaboration prefers the 
use of cross-spectra even at a single frequency to gain ad¬ 
ditional immunity to systematics and to avoid the need 
to noise debias auto-spectra. The TT spectrum is higher 
than ACDM around £ = 200 -presumably due to a dust 
contribution. The EE and BB spectra are noisy, but 
both appear to show an excess over ACDM for £ < 150— 
again presumably due to dust. We note that these spec¬ 
tra do not appear to follow the power-law expectation 
mentioned in Sec. IIB, but we emphasize that the error 
bars contain no sample variance on any dust component 
(Gaussian or otherwise). 

The center column of Fig. 2 shows cross-spectra be¬ 
tween BICEP2/KecA; and Planck maps. For TE one 
can use the T-modes from BICEP2 and the E-modes 
from Planck or vice versa and both options are shown. 
Since the T-modes are very similar between the two ex¬ 
periments, these TE spectra look similar to the single¬ 
experiment TE spectrum which shares the E-modes. 
The EE and BB cross-spectra are the most interesting— 
there appears to be a highly significant detection of cor¬ 
related E-mode power between 150 and 353 GHz, with 
the pattern being much brighter at 353, consistent with 
the expectation from dust. We also see hints of detection 
in the EE spectrum- -while dust E-modes are subdomi¬ 
nant to the cosmological signal at 150 GHz, the weak dust 
contribution enhances the BK150xP353 cross-spectrum 
at £ « 100. 

The polarized dust SED model mentioned in Sec. IIB 
implies that dust emission is approximately 25 times 
brighter in the Planck 353 GHz band than it is in the 
BICEP2/EecA; 150 GHz band (integrating appropriately 
over the instrumental bandpasses). The expectation for a 
dust-dominated spectrum is thus that the BK150xP353 
cross-spectrum should have an amplitude 25 times that 
of BK150xBK150, and P353xP353 should be 25 times 
higher again. The y-axis scaling in the bottom row of 
Fig. 2 has been adjusted so that a dust signal obeying this 
rule will have equal apparent amplitude in each panel. 
We see that a substantial amount of the BK150xBK150 
signal indeed appears to be due to dust. 

To make a rough estimate of the significance of devi¬ 
ation from lensed-ACDM, we calculate y 2 and \ (sum 
of deviations) for each of the EE and BB spectra and 
show these in Fig. 2. For the nine bandpowers used the 
expectation value/standard-deviation for % 2 and y are 
9/4.2 and 0/3 respectively. We see that BK150xBK150 
and BK150xP353 are highly significant in EE, while 
P353xP353 has modest significance in both EE and EE. 


Figure 3 shows EE and EE cross-spectra between BI- 
CEP2 /Keck and all of the polarized frequencies of Planck 
(also including the BICEP2/Kec£' auto-spectra). For the 
five bandpowers shown the expectation value/standard- 
deviation for x 2 and X are 5/3.1 and 0/2.2 respectively. 
As already noted, the BK150xBK150 and BK150xP353 
EE spectra show highly significant excesses. Addition¬ 
ally, there is evidence for excess EE in BK150xP217 
spectrum, and for excess EE in BK150xP353. The other 
spectra in Fig. 3 show no strong evidence for excess, al¬ 
though we note that only one of the \ values is negative. 
There is weak evidence for excess in the BK150xP70 EE 
spectrum but none in BK150xP30 so this is presumably 
just a noise fluctation. 

There are a large number of additional Planck- only 
spectra, which are not plotted here. The noise on these is 
large and all are consistent with ACDM, with the possible 
exception of P217xP353, where modest evidence for an 
excess is seen in both EE and EE (see e.g., Figure 10 of 
PIP-XXX). 

D. Consistency of BICEP2 and Keck Array spectra 

The EE auto-spectra for BICEP2 and Keck Array in 
the lower left panel of Fig. 2 appear to differ by more 
than might be expected, given that the BICEP2 and 
Keck maps cover almost exactly the same region of sky. 
However, the error bars in this figure are the standard 
deviations of lensed-ACDM+noise simulations; while the 
signal is largely common between the two experiments 
the noise is not, and the signal-noise cross terms produce 
substantial additional fluctuation of the difference. The 
correct way to quantify this is to compare the difference 
of the real data to the pair-wise differences of simulations, 
using common input skies that have power similar to that 
observed in the real data. This was done in Section 8 of 
BK-V and the BICEP2 and Keck maps were shown to be 
statistically compatible. In an analogous manner we can 
also ask if the B150xP353 and K150xP353 EE cross¬ 
spectra shown in the bottom middle panel of Fig. 2 are 
compatible. Figure 4 shows the results. We calculate the 
X 2 a nd y statistics on these difference spectra and com¬ 
pare to the simulated distributions exactly as in BK-V. 
The probability to exceed (PTE) the observed values is 
given in the figure for bandpowers 1-5 (20 < l < 200) 
and 1 9 (20 < i < 330). There is no evidence that these 
spectra are statistically incompatible. 

E. Alternative power spectrum estimation 

We check the reliability of the power spectrum estima¬ 
tion with an alternative pipeline. The filtered and puri¬ 
fied Planck and BICEP2 maps used to make the spectra 
shown in Fig. 2 are transformed back into the HEALPix 
pixelization using cubic spline interpolation. The E- 
mode cross-power is then computed with the Xpol [39] 
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FIG. 2. Single- and cross-frequency spectra between BICEP2 /Keck maps at 150 GHz and Planck maps at 353 GHz. The red 
curves show the lensed-ACDM expectations. The left column shows single-frequency spectra of the BICEP2, Keck Array and 
combined BICEP2/Aecfc maps. The BICEP2 spectra are identical to those in BK-I, while the Keck Array and combined are 
as given in BK-V. The center column shows cross-frequency spectra between BICEP2/ Keck maps and Planck 353 GHz maps. 
The right column shows Planck 353 GHz data-split cross-spectra. In all cases the error bars are the standard deviations of 
lensed-ACDM+noise simulations and hence contain no sample variance on any other component. For EE and BB the y 2 and 
X (sum of deviations) versus lensed-ACDM for the nine bandpowers shown is marked at upper/lower left (for the combined 
BICEP2 /Keck points and DSlxDS2). In the bottom row (for BB) the center and right panels have a scaling applied such 
that signal from dust with the fiducial frequency spectrum would produce signal with the same apparent amplitude as in the 
150 GHz panel on the left (as indicated by the right-side y- axes). We see from the significant excess apparent in the bottom 
center panel that a substantial amount of the signal detected at 150 GHz by BICEP2 and Keck Array indeed appears to be due 
to dust. 


and PtjreCl [40] estimators. Figure 5 shows the differ¬ 
ence between these alternative bandpowers and the stan¬ 
dard bandpowers for the B150xP353 BB cross spectrum. 
As in Fig. 4 the errorbars are the standard deviations of 
pairwise differences of simulations, which share common 
input skies and have power similar to that observed in 
the real data. The agreement is not expected to be ex¬ 
act due to the differing bandpower window functions, but 
the differences of the real bandpowers are consistent with 
those of the simulations. 


III. LIKELIHOOD ANALYSIS 
A. Algorithm 

While it is conventional in plots like Fig. 2 to present 
bandpowers with symmetric error bars, it is important to 
appreciate that this is an approximation. The likelihood 
of an observed bandpower for a given model expecta¬ 
tion value is generally an asymmetric function, which can 
be computed given knowledge of the noise level(s). To 
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FIG. 4. Differences of B150xP353 and K150xP353 BB cross¬ 
spectra. The error bars are the standard deviations of the 
pairwise differences of signal+noise simulations that share 
common input skies. The probability to exceed the observed 
values of y 2 an d \ statistics, as evaluated against the simula¬ 
tions, is quoted for bandpower ranges 1-5 (20 < l < 200) and 
1-9 (20 < l < 330). There is no evidence that these spectra 
are statistically incompatible. 


• B150xP353-K150xP353 


PTEs: 1 -5 X 2 0.273, 1 -9 X 2 0.415, 1 -5 X 0.802, 1 -9 X 0.425 
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FIG. 5. Differences of B150xP353 BB cross-spectra from the 
standard power spectrum estimator and alternate estimators. 
The error bars are the standard deviations of the pairwise dif¬ 
ferences of signal+noise simulations that share common input 
skies. The probability to exceed the observed values of \ 2 and 
X statistics, as evaluated against the simulations, is quoted for 
bandpower ranges 1 5 (20 < £ < 200) and 1 -9 (20 < £ < 330). 
We see that the differences of the real spectra are consistent 
with the differences of the simulations. 
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FIG. 3. EE (left column) and BB (right column) cross¬ 
spectra between BICEP2/A'ecfc maps and all of the polar¬ 
ized frequencies of Planck. In all cases the quantity plotted is 
£(£+l)Ci/2n in units of +Kq MB , and the red curves show the 
lensed-ACDM expectations. The error bars are the standard 
deviations of lensed-ACDM+noise simulations and hence con¬ 
tain no sample variance on any other component. Also note 
that the y -axis scales differ from panel to panel in the right 
column. The y 2 an d X (sum of deviations) versus lensed- 
ACDM for the five bandpowers shown is marked at upper left. 
There are no additional strong detections of deviation from 
lensed-ACDM over those already shown in Fig. 2 although 
BK150xP217 shows some evidence of excess. 


compute the joint likelihood of an ensemble of measured 
bandpower values it is of course necessary to consider 
their full covariance—this is especially important when 
using spectra taken at different frequencies on the same 
field, where the signal covariance can be very strong. 

We compute the bandpower covariance using full 
simulations of signal-cross-signal, noise-cross-noise, and 
signal-cross-noise. From these, we can construct the co- 
variance matrix for a general model containing multiple 






































































































































10 


signal components with any desired set of SEDs. When 
we do this we deliberately exclude terms whose expecta¬ 
tion value is zero, in order to reduce noise in the resulting 
matrix due to the limited number of simulated realiza¬ 
tions. 

To compute the joint likelihood of the data for any 
given proposed model we use the Hamimeche-Lewis [41] 
approximation (HL; see Section 9.1 of Ref. [42] for imple¬ 
mentation details). Here we extend the method to deal 
with single- and cross-frequency spectra, and the covari¬ 
ances thereof, in an analogous manner to the treatment 
of, for example, TT, TE, and EE in the standard HL 
method. The HL formulation requires that the band- 
power covariance matrix be determined for only a single 
“fiducial model.” We compute multi-dimensional grids 
of models explicitly and/or use COSMOMC [43] to sample 
the parameter space. 


B. Fiducial analysis 

As an extension of the simplest lensed-ACDM 
paradigm, we initially consider a two component model 
of IGW with amplitude r, plus dust with amplitude A<j 
(specified at 353 GHz and t = 80). (Here we assume 
that the spectral index of the tensor modes (rit) is zero, 
and a scalar pivot scale of 0.05 Mpc -1 ; all values of r 
quoted in this paper are ro.os unless noted otherwise.) 
Figure 6 shows the results of fitting such a model to 
BB bandpowers taken between BICEP2 /Keck and the 
217 and 353 GHz bands of Planck , using bandpowers 1-5 
(20 < £ < 200). For the Planck single-frequency case, 
the cross-spectrum of detector-sets (DSlxDS2) is used, 
following PIP-XXX. The dust is modeled as a power law 
T>t oc £ -0 - 42 , with free amplitude Ad and scaling with 
frequency according to the modified blackbody model. 

As discussed in Sec. IIB the simple modified black- 
body model is shown empirically in Ref. [26] to describe 
the mean polarized dust SED at mid Galactic latitudes 
to an accuracy of a few percent over the frequency range 
100-353 GHz, with variation of the /3d parameter being 
sufficient to characterize the patch-to-patch variation. 
Since it is not possible to constrain /3d using the BI- 
CEP2 /Keck/ and Planck cross spectral bandpowers alone 
a tight Gaussian prior /3d = 1.59±0.11 is imposed, the un¬ 
certainty being scaled from the observed patch-to-patch 
variation at intermediate Galactic latitudes in Ref. [26], 
as explained in PIP-XXX. This prior assumes that the 
SED of dust polarization at intermediate latitudes [26] 
applies to the high latitude BICEP2 /Keck field. From 
dust astrophysics, we expect variations of the dust SED 
in intensity and polarization to be correlated [18]. We 
thus tested our assumption by measuring the /3d of the 
dust total intensity in the BICEP2 /Keck field using the 
template fitting analysis described in Ref. [44], and find 
the same value. 

In Fig. 6 we see that the BICEP2 data produce an r 
likelihood that peaks higher than that for the Keck Ar¬ 


ray data. This is because for £ < 120 the auto-spectrum 
B150xB150 is higher than for K150xK150, while the 
cross-spectrum B150xP353 is lower than K150xP353 
(see Fig. 2). However, recall that both pairs of spectra 
B150xB150/K150xK150 and B150xP353/K150xP353 
have been shown to be consistent within noise fluctua¬ 
tion (see Sec. IID). In Sec. IV A these likelihood results 
are also found to be compatible. Given the consistency 
between the two experiments, the combined result gives 
the best available measurement of the sky. 

The combined curves (BK+P) in the left and cen¬ 
ter panels of Fig. 6 yield the following results: r = 
0.048j/p;Q32, r < 0.12 at 95% confidence, and Ad = 
3.3+0 g. For r the zero-to-peak likelihood ratio is 0.38. 
Taking \(l - f (-2 log L 0 /L peak )), where / is the y 2 cdf 
(for one degree of freedom), we estimate that the proba¬ 
bility to get a number smaller than this is 8% if in fact 
r = 0. For Ad the zero-to-peak ratio is 1.8 x 10~ 6 corre¬ 
sponding to a smaller-than probability of 1.4 x 10 -7 , and 
a 5.1(J detection of dust power. 

The maximum likelihood model on the grid has pa¬ 
rameters r = 0.05, Ad = 3.30 /xK 2 (and /3d = 1.6). Com¬ 
puting the bandpower covariance matrix for this model, 
we obtain a y 2 of 40.8. Using 28 degrees of freedom— 
5 bandpowers times 6 spectra, minus 2 fit parameters 
(since /3d is not really free)—gives a PTE of 0.06. The 
largest contributions to y 2 come from the P353xP353 
spectrum shown in the lower right panel of Fig. 2. 

C. Variations from the fiducial data set and model 

We now investigate a number of variations from the 
fiducial analysis to see what difference these make to the 
constraint on r. 

• Choice of Planck single-frequency spectra: 

switching the Planck single-frequency spectra to 
use one of the alternative data splits (yearly or half¬ 
ring instead of detector-set) makes little difference 
(see Fig. 7). 

• Using only 150 and 353 GHz: dropping the 
spectra involving 217 GHz from consideration also 
has little effect (see Fig. 7). 

• Using only BK150xBK150 and 

BK150xP353: also excluding the 353 GHz 

single-frequency spectrum from consideration 
makes little difference. The statistical weight of 
the BK150xBK150 and BK150xP353 spectra 
dominate (see Fig. 7). 

• Extending the bandpower range: going back 
to the base data set and extending the range of 
bandpowers considered to 1-9 (corresponding to 
20 < t < 330) makes very little difference—the 
dominant statistical weight is with the lower band- 
powers (see Fig. 7). 
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FIG. 6. Likelihood results from a basic lensed-ACDM+r+dust model, fitting BB auto- and cross-spectra taken between maps 
at 150 GHz, 217, and 353 GHz. The 217 and 353 GHz maps come from Planck. The primary results (heavy black) use the 
150 GHz combined maps from BICEP2 /Keck. Alternate curves (light blue and red) show how the results vary when the 
BICEP2 and Keck Array only maps are used. In all cases a Gaussian prior is placed on the dust frequency spectrum parameter 
/3a = 1.59 ± 0.11. In the right panel the two dimensional contours enclose 68% and 95% of the total likelihood. 


• Including EE spectra: we can also include in 
the fits the EE spectra shown in Fig. 3. PIP-XXX 
(figures 5 and A.3) shows that the level of EE from 
Galactic dust is on average around twice the level 
of BB. However, there are substantial variations 
in this ratio from sky-patch to sky-patch. Setting 
EE/BB = 2 we find that the constraint on Ad 
narrows, while the r constraint changes little; this 
latter result is also shown in Fig. 7. The maximum 
likelihood model on the grid is unchanged and its 
X 2 PTE is acceptable. 

• Relaxing the /3d prior: relaxing the prior on the 
dust spectral index to fid = 1-59 ± 0.33 pushes the 
peak of the r constraint up (see Fig. 7). However, 
it is not clear if this looser prior is self consistent; 
if the frequency spectral index varied significantly 
across the sky it would invalidate cross-spectral 
analysis, but there is strong evidence against such 
variation at high latitude, as explained in Sec. V A. 
Nevertheless, it is important to appreciate that the 
r constraint curves shown in Fig. 6 shift left (right) 
when assuming a lower (higher) value of fid- For 
/3d = 1-3 ± 0.11 the peak is at r = 0.021 and for 
/3d = 1-9 ± 0.11 the peak is at r = 0.073. 

• Varying the dust power spectrum shape: in 

the fiducial analysis the dust spatial power spec¬ 
trum is assumed to be a power law with Ve oc 
t~ 0A2 . Marginalizing over spectral indices in the 
range —0.8 to 0 we find little change in the r con¬ 
straint (see also Sec. IV B for an alternate relax¬ 
ation of the assumptions regarding the spatial prop¬ 
erties of the dust pattern). 

• Using Gaussian determinant likelihood: the 

fiducial analysis uses the HL likelihood approx¬ 
imation, as described in Sec. Ill A. An alterna¬ 
tive is to recompute the covariance matrix C at 
each point in parameter space and take L = 


det (C) -1 / 2 exp (— (d T C _1 d)/2), where d is the de¬ 
viation of the observed bandpowers from the model 
expectation values. This results in an r constraint 
which peaks slightly lower, as shown in Fig. 7. Run¬ 
ning both methods on the simulated realizations 
described in Sec. IV A, indicates that such a dif¬ 
ference is not unexpected and that there may be 
a small systematic downward bias in the Gaussian 
determinant method. 

• Varying the HL fiducial model: as mentioned 
in Sec. Ill A the HL likelihood formulation re¬ 
quires that the expectation values and bandpower 
covariance matrix be provided for a single “fidu¬ 
cial model” (not to be confused with the “fidu¬ 
cial analysis” of Sec. IIIB). The results from HL 
are supposed to be rather insensitive to the choice 
of this model, although preferably it should be 
close to reality. Normally we use the lensed- 
ACDM+dust simulations described in Sec. IV A be¬ 
low. Switching this to lensed-ACDM+r=0.2 pro¬ 
duces no change on average in the simulations, al¬ 
though it does cause any given realization to shift 
slightly—the change for the real data case is shown 
in Fig. 7. 

• Adding synchrotron: BK-I took the WMAP K- 
band (23 GHz) map, extrapolated it to 150 GHz ac¬ 
cording to v~ 33 (mean value within the BICEP2 
field of the MCMC “Model f” spectral index map 
provided by WMAP [2]), and found a negligible 
predicted contribution (r synCj i 50 = 0.00 08±0.0041). 
Figure 3 does not offer strong motivation to reex¬ 
amine this finding—the only significant detections 
of correlated BB power are in the BK150xP353 
and, to a lesser extent, BK150xP217 spectra. How¬ 
ever, here we proceed to a fit including all the 
polarized bands of Planck (as shown in Fig. 3) 
and adding a synchrotron component to the base 
lensed-ACDM+noise+r+dust model. We take syn- 
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chrotron to have a power law spectrum T>i oc 
£- 0.6 [ 23 ] ; with free amplitude A sync , where A sync is 
the amplitude at l = 80 and at 150 GHz, and scal¬ 
ing with frequency according to v~ 33 . In such a 
scenario we can vary the degree of correlation that 
is assumed between the dust and synchrotron sky 
patterns. Figure 8 shows results for the uncorre¬ 
lated and fully correlated cases. Marginalizing over 
r and Ad we find A sync < 0.0003 ^K 2 at 95% con¬ 
fidence for the uncorrelated case, and many times 
smaller for the correlated. This last is because once 
one has a detection of dust it effectively becomes 
a template for the synchrotron. This synchrotron 
limit is driven by the Planck 30 GHz band -we ob¬ 
tain almost identical results when adding only this 
band, and a much softer limit when not including it. 
If we instead assume synchrotron scaling of i/~ 3 0 
the limit on A sync is approximately doubled for the 
uncorrelated case and reduced for the correlated. 
(Because the DSlxDS2 data-split is not available 
for the Planck LFI bands we switch to YlxY2 for 
this variant analysis, and so we compare to this 
case in Fig. 8 rather than the usual fiducial case.) 

• Varying lensing amplitude: in the fiducial anal¬ 
ysis the amplitude of the lensing effect is held 
fixed at the ACDM expectation (Al = 1). Using 
their own and other data, the Planck Collabora¬ 
tion quote a limit on the amplitude of the lens¬ 
ing effect versus the ACDM expectation of A L = 
0.99 ± 0.05 [3]. Allowing Al to float freely, and 
using all nine bandpowers, we obtain the results 
shown in Fig. 9—there is only weak degeneracy be¬ 
tween Al and both r and A j. Marginalizing over r 
and Ad we find Al = 1.13 ± 0.18 with a likelihood 
ratio between zero and peak of 3 x 10 -11 . Using the 
expression given in Sec. IIIB this corresponds to a 
smaller-than probability of 2 x 10~ 12 , equivalent to 
a 7.0 a detection of lensing in the BB spectrum. 
We note this is the most significant to-date direct 
measurement of lensing in H-rnode polarization. 


IV. LIKELIHOOD VALIDATION 

A. Validation with simulations 

We run the algorithm used in Sec. Ill B on ensembles of 
simulated realizations to check its performance. We first 
consider a model where r = 0 and Ad = 3.6 /iK 2 , this lat¬ 
ter being close to the value favored by the data in a dust- 
only scenario [45]. We generate Gaussian random real¬ 
izations using the fiducial spatial power law T>e oc £~ 0 42 , 
scale these to the various frequency bands using the mod¬ 
ified blackbody law with Td = 19.6 K and /3d = 1-59, 
and add to the usual realizations of lensed-ACDM+noise. 
Figure 10 shows some of the resulting r and Ad constraint 



Fiducial analysis 
Y1xY2 
no 217GHz 

Only BKxBK&BKxP353 
9 bandpowers 
Inc. EE (EE/BB=2) 
relax P d prior 

Gauss det 
alt. HL fid. model 


FIG. 7. Likelihood results when varying the data sets used 
and the model priors—see Sec III C for details. 



FIG. 8. Likelihood results for a fit when adding the lower 
frequency bands of Planck, and extending the model to in¬ 
clude a synchrotron component. The results for two differ¬ 
ent assumed degrees of correlation between the dust and syn¬ 
chrotron sky patterns are compared to those for the compa¬ 
rable model without synchrotron (see text for details). 


curves, with the result for the real data from Fig. 6 over¬ 
plotted. As expected, approximately 50% of the r like¬ 
lihoods peak above zero. The median 95% upper limit 
is r < 0.075. We find that 8% of the realizations have 
a ratio Lo/^peak less than the 0.38 observed in the real 
data, in agreement with the estimate in Sec. IIIB. Run¬ 
ning these dust-only realizations for BICEP2 only and 
Keck Array only, we find that the shift in the maximum 
likelihood value of r seen in the real data in Fig. 6 is 
exceeded in about 10% of the simulations. 

The above simulations assume that the dust compo¬ 
nent follows on average the fiducial T>t oc £~ 0A2 spatial 
power law, and fluctuates around it in a Gaussian man¬ 
ner. To obtain sample dust sky patterns that may deviate 
from this behavior in a way which better reflects reality, 
we take the pre-launch version of the Planck Sky Model 
(PSM; version 1.7.8 run in “simulation” mode) [24] eval¬ 
uated in the Planck 353 GHz band and pull out the same 
352 \b\ >35° partially overlapping regions used in PIP- 
XXX. We then scale these to the other bands and proceed 
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FIG. 9. Likelihood results for a fit allowing the lensing 
scale factor A l to float freely and using all nine bandpowers. 
Marginalizing over r and Ad, we find that Al = 1.13 ± 0.18 
and A l = 0 is ruled out with 7.0 a significance. 




FIG. 10. Likelihoods for r and Ad, using BICEP2 /Keck 
and Planck, as plotted in Fig. 6, overplotted on constraints 
obtained from realizations of a lensed-ACDM+noise+dust 
model with dust power similar to that favored by the real 
data (Ad = 3.6/rK 2 ). Half of the r curves peak at zero as 
expected. 


as before. Some of the regions have dust power orders of 
magnitude higher than the real data and we cut them 
out (selecting 139 regions with peak Ad < 20/uK 2 ). Fig¬ 
ure 11 presents the results. The r likelihoods will broaden 
as the level of Ad increases, and we should therefore not 
be surprised if the fraction of realizations peaking at a 
value higher than the real data is increased compared to 
the simulations with mean Ad = 3.6/rK 2 . However we 
still expect that on average 50% will peak above zero and 
approximately 8% will have an Lo/L pea k ratio less than 
the 0.38 observed in the real data. In fact we find 57% 
and 7%, respectively, consistent with the expected val¬ 
ues. There is one realization which has a nominal (false) 
detection of non-zero r of 3.3<r, although this turns out to 
also have one of the lowest L 0 /L pea k ratios in the Gaus¬ 
sian simulations shown in Fig. 10 (with which it shares 
the CMB and noise components), so this is apparently 
just a relatively unlikely fluctuation. 
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FIG. 11. Constraints obtained when adding dust realizations 
from the Planck Sky Model version 1.7.8 to the base lensed- 
ACDM+noise simulations. (Curves for 139 regions with peak 
Ad < 20 fiK 2 are plotted.) We see that the results for r 
are unbiased in the presence of dust realizations which do 
not necessarily follow the £~ 0 - 42 power law or have Gaussian 
fluctuations about it. 


B. Subtraction of scaled spectra 

As previously mentioned, the modified blackbody 
model predicts that dust emission is 4% as bright in the 
BICEP2 band as it is in the Planck 353 GHz band. There¬ 
fore, taking the auto- and cross-spectra of the combined 
BICEP2/Aecfc maps and the Planck 353 GHz maps, as 
shown in the bottom row of Fig. 2, and evaluating 
(BK x BK — aBK x P)/(l — a), at a — agd cleans out the 
dust contribution (where ckfid = 0.04). The upper panel 
of Fig. 12 shows the result. 

As an alternative to the full likelihood analysis pre¬ 
sented in Sec. IIIB, we can instead work with the dif¬ 
ferenced spectra from above, a method we denote the 
“cleaning” approach. If agd were the true value, the ex¬ 
pectation value of this combination over CMB and noise 
would have no dust contribution. However, dust would 
still contribute to its variance, but only through its 2- 
point function. In practice, we do not know a perfectly, 
and this uncertainty needs to be accounted for in a like¬ 
lihood constructed from the differenced spectra. Our ap¬ 
proach is to treat the differenced spectra as a form of 
data compression, and to compute the expectation value 
as a function of r, Ad, and /3d at each point in parameter 
space (the dust dependence enters for a (/3d) 7 ^ afid)- We 
use the method of Ref. [41], with a fiducial covariance 
matrix, to build a likelihood for the difference spectra, 
and marginalize over Ad and /3d, and hence a, adopting 
the prior /3d = 1.59±0.11. This alternative likelihood has 
the advantage of being less sensitive to non-Gaussianity 
of the dust, since only the 2-point function of the dust 
affects the covariance of the differenced spectra close to 
add) while the full analysis may, in principle, be affected 
by the non-Gaussianity of the dust through 4-point con¬ 
tributions to power spectra covariances. This cleaning 
approach does, however, ignore the (small amount of) ad¬ 
ditional information available at other frequencies. The 
lower panel of Fig. 12 compares the result to the fiducial 


















14 




r 


FIG. 12. Upper: BB spectrum of the BICEP2/A'ecfcmaps be¬ 
fore and after subtraction of the dust contribution, estimated 
from the cross-spectrum with Planck 353 GHz. The error bars 
are the standard deviations of simulations, which, in the lat¬ 
ter case, have been scaled and combined in the same way. The 
inner error bars are from lensed-ACDM+noise simulations as 
in the previous plots, while the outer error bars are from 
the lensed-ACDM+noise+dust simulations. The red curve 
shows the lensed-ACDM expectation. Lower: constraint on r 
derived from the cleaned spectrum compared to the fiducial 
analysis shown in Fig. 6. 


analysis with the full multi-spectra likelihood. It is clear 
from the widths of the likelihood curves that compressing 
the spectra to form the cleaned difference results in very 
little loss of information on r. The difference in peak 
values arises from the different data treatments and is 
consistent with the scatter seen across simulations. Fi¬ 
nally, we note that one could also form a combination 
(BK x BK — 2aBK x P + a 2 P x P)/(l — a) 2 in which dust 
does not enter at all for a = agd- However, the variance 
of this combination of spectra is large due to the Planck 
noise levels, and likelihoods built from this combination 
are considerably less constraining. 


V. POSSIBLE CAUSES OF DECORRELATION 

Any systematic error that suppresses the BK150xP353 
cross-frequency spectrum with respect to the 
BK150xBK150 and P353xP353 single-frequency 
spectra would cause a systematic upward bias on the r 
constraint. Here we investigate a couple of possibilities. 

A. Spatially varying dust frequency spectrum 

If the frequency dependence of polarized dust emission 
varied from place to place on the sky, it would cause the 
150 GHz and 353 GHz dust sky patterns to decorrelate 
and suppress the BK150xP353 cross-frequency spectrum 
relative to the single-frequency spectra. The assump¬ 
tion made so far in this paper is that such decorrela¬ 
tion is negligible. In fact PIP-XXX implicitly tests for 
such variation in their Figure 6, where the Planck single- 
and cross-frequency spectra are compared to the modi¬ 
fied blackbody model (with the cross-frequency spectra 
plotted at the geometric mean of their respective frequen¬ 
cies). This plot is for an average over a large region of low 
foreground sky (24%); however, note that if there were 
spatial variation of the spectral behavior anywhere in this 
region it would cause suppression of the cross-frequency 
spectra with respect to the single-frequency spectra. 

PIP-XXX also tests explicitly for evidence of decorre¬ 
lation of the dust pattern across frequencies. Their fig¬ 
ure E.l shows the results for large and small sky patches. 
The signal-to-noise ratio is low in clean regions, but no 
evidence of decorrelation is found. 

As a further check, we artificially suppress the ampli¬ 
tude of the BK150 xP353 spectra in the Gaussian dust- 
only simulations (see Sec. IV A) by a conservative 10% 
(PIP-XXX sets a 7% upper limit). We find that the 
maximum likelihood value for r shifts up by an average 
of 0.018, while Ad shifts down by an average of 0.43 pK 2 , 
with the size of the shift proportional to the magnitude of 
the dust power in each given realization. This behavior 
is readily understandable—since the BK150xBK150 and 
BK150xP353 spectra dominate the statistical weight, a 
decrease of the latter is interpreted as a reduction in dust 
power, which is compensated by an increase in r. The 
bias on r will be linearly related to the assumed decorre¬ 
lation factor. 


B. Calibration, analysis etc. 

Figure 3 shows that the EE spectrum BK150xBK150 
is extremely similar to that for BK150xP143. We 
can compare such spectra to set limits on possible 
decorrelation between the BICEP2/AVcfc and Planck 
maps arising from any instrumental or analysis re¬ 
lated effect, including differential pointing, polarization 
angle mis-characterization, etc. Taking the ratio of 
BK150xP143 to the geometric mean of BK150xBK150 
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and P143HlxP143H2, we find that for TT the decorre¬ 
lation is approximately 0.1%. For EE the signal-to-noise 
ratio is lower, but decorrelation is limited to below 2%, 
and consistent with zero when compared to the fluctua¬ 
tion of signal+noise simulations. 


VI. CONCLUSIONS 

BK-I reported a highly significant detection of 13-mode 
polarization, at 150 GHz, in excess of the lensed-ACDM 
expectation over the range 30 < £ < 150. This excess has 
been confirmed by additional data on the same field from 
the sucessor experiment Keck Array. PIP-XXX found 
that the level of dust power in a field centered on the 
PAGEP2 / Keck region (but somewhat larger than it) is of 
the same magnitude as the reported excess, but noted 
that, “the present uncertainties are large,” and that a 
joint analysis was required. 

In this paper we have performed this joint analysis, us¬ 
ing the combined PAGEP2 / Keck maps. Cross-correlating 
these maps against all of the polarized frequency bands 
of Planck we find a highly significant 13-mode detection 
only in the cross spectrum with 353 GHz. We emphasize 
that this 150 x 353 GHz cross-spectrum has a much higher 
signal-to-noise ratio than the 353 GHz single-frequency 
spectrum that PIP-XXX analyzed. 

We have analyzed the data using a multi-frequency, 
multi-component fit. In this fit it is necessary to im¬ 
pose a prior on the variation of the brightness of the 
polarized dust emission with observing frequency, since 
the available data are unable to constrain this alone, due 
to the relatively low signal-to-noise ratio in 13-mode po¬ 
larization at 353 GHz. However based on the available 
information from Planck on the frequency dependence of 
polarized dust emission across the mid and high Galac¬ 
tic latitude sky, and the patch-to-patclr stability thereof, 
this prior appears to be justified and conservative. 

We have shown that the final constraint on the tensor- 
to-scalar ratio r is very stable when varying the frequency 
bands used, as well as the model priors. The result does 
differ when using the BICEP2 and Keck Array data alone 
rather than in combination, but the difference is com¬ 
patible with noise fluctuation. Expanding the model to 
include synchrotron emission, while also including lower 
Planck frequencies, does not change the result. 

Allowing the amplitude of lensing to be free, we ob¬ 
tain Al = 1.13 ± 0.18, with a significance of detection of 
7.Ocr. This is the most significant direct detection to-date 
of lensing in i3-mode polarization, even compared to ex¬ 
periments with higher angular resolution. The POLAR- 
BEAR experiment has reported a detection of 33-mode 
lensing on smaller angular scales (500 < t < 2100), re¬ 
jecting the Al = 0 hypothesis at 97.2% confidence [46]. 
Additionally, ACT [47] and SPT [48] have reported lens¬ 
ing detections in polarization in cross-correlation with 
some other tracer of the dark matter distribution on the 
sky. 


We have validated the main likelihood analysis on sim¬ 
ulations of a dust-only model and performed a simple 
subtraction of scaled spectra, which approximates a map- 
based dust cleaning (obtaining an r constraint curve that 
peaks somewhat lower). Finally we investigated the pos¬ 
sibility of astrophysical or instrumental decorrelation of 
the sky patterns between experiments or frequencies and 
find no evidence for relevant bias. 

The final result is expressed as a likelihood curve for 
r, and yields an upper limit r < 0.12 at 95% confidence. 
The median limit in the lensed-ACDM+noise+dust sim¬ 
ulations is r < 0.075. It is interesting to compare this 
latter to dust-free simulations using only BICEP2/Aecfc 
where the median limit is r < 0.03—the difference rep¬ 
resents the limitation due to noise in the Planck maps, 
when marginalizing over dust. The r constraint curve 
peaks at r = 0.05 but disfavors zero only by a fac¬ 
tor of 2.5. This is expected by chance 8% of the time, 
as confirmed in simulations of a dust-only model. We 
emphasize that this significance is too low to be inter¬ 
preted as a detection of primordial 33-modes. Transform¬ 
ing the Planck temperature-only 95% confidence limit of 
f 0.002 < 0.11 [3] to the pivot scale used in this paper 
yields ro.os < 0.12, compatible with the present result. 

A C0SM0MC module containing the bandpowers for all 
cross-spectra between the combined BICEP2/Aecfc maps 
and all of the frequencies of Planck is available for down¬ 
load at http://bicepkeck.org. 

In order to further constrain or detect IGW, addi¬ 
tional data are required. The Planck Collaboration may 
be able to make progress alone using the large angular 
scale “reionization bump,” if systematics can be appro¬ 
priately controlled [49]. To take small patch “recombi¬ 
nation bump” studies of the type pursued here to the 
next level, data with signal-to-noise comparable to that 
achieved by BICEP2/Aecfc at 150 GHz are required at 
more than one frequency. Figure 13 summarizes the sit¬ 
uation. The BICEP2 /Keck noise is much lower in the BI- 
CEiP2jKeck field than the Planck noise. However, since 
dust emission is dramatically brighter at 353 GHz, it is 
detected in the cross-spectrum between BICEP2 /Keck 
and Planck 353 GHz. Synchrotron is not detected and 
the crossover frequency with dust is < 100 GHz. Planck’s 
PR2 data release [50] shows that for the cleanest 73% of 
the sky, at 40’ scales, the polarized foreground minimum 
is at ~80-90 GHz. During the 2014 season, two of the 
Keck Array receivers observed in the 95 GHz band and 
these data are under active analysis. BICEP3 will add 
substantial additional sensitivity at 95 GHz in the 2015, 
and especially 2016, seasons. Meanwhile many other 
ground-based and sub-orbital experiments are making 
measurements at a variety of frequencies and sky cov¬ 
erage fractions. 
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Nominal band center [GHz] 


FIG. 13. Expectation values, and uncertainties thereon, 
for the i ~ 80 BB bandpower in the BICEP2 /Keck field. 
The green and magenta lines correspond to the expected sig¬ 
nal power of lensed-ACDM and r = 0.05. Since CMB units 
are used, the levels corresponding to these are flat with fre¬ 
quency. The grey band shows the best fit dust model (see Sec¬ 
tion III B) and the blue shaded region shows the allowed re¬ 
gion for synchrotron (see Sec. IIIC). The BICEP2/Aecfc noise 
uncertainty is shown as a single starred point, and the noise 
uncertainties of the Planck single-frequency spectra evaluated 
in the BICEP2 /Keck held are shown in red. The blue points 
show the noise uncertainty of the cross-spectra taken between 
BICEP2/Kecfc and, from left to right, Planck 30, 44, 70, 100, 
143, 217 & 353 GHz, and plotted at horizontal positions such 
that they can be compared vertically with the dust and sync 
curves. 
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